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Q I A new type of object called "Very Low Luminosity Objects (VeLLOs)" has 

^ I been discovered by the Spitzer Space Telescope. VeLLOs might be substellar 
objects forming by accretion. However, some VeLLOs are associated with strong 

r-| I outflows, indicating the previous existence of massive accretion. The thermal 

^1 history, which significantly affects the chemistry, between substellar objects with 

O I a continuous low accretion rate and objects in a quiescent phase after massive ac- 

^ [ cretion (outburst) must be greatly different. In this study, the chemical evolution 
has been calculated in an episodic accretion model to show that CO and N2H"'" 

_ have a relation different from starless cores or Class O/I objects. Furthermore, 

> I the CO2 ice feature at 15.2 /xmm will be a good tracer of the thermal process in 

^ ■ VeLLOs. 
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INTRODUCTION 



j_| ' The Spitzer Space Telescope (SST), with its very high mid-infrared sensitivity, has re- 

^ vealed a new type of low-mass protostars called "Very Low Luminosity Objects (VeLLOs; 

e.g.. Young et al. 2004). Defining the internal luminosity of a source. Lint, to be the total 
luminosity of the central protostar and circumstellar disk (if present), a VeLLO is defined 
to be an object embedded within a dense core that meets the criterion Li^t < 0.1 Lq (Di 
Francesco et al. 2007). VeLLOs are primarily being discovered in cores previously classified 
as starless based on observations with the Infrared Astronomical Satellite (IRAS) and the 
Infrared Space Observatory (ISO). The three most well-studied VeLLOs to date are IRAM 
04191 (Andre et al. 1999; Dunham et al. 2006), L1014 (Young et al. 2004; Crapsi et al. 
2005; Bourke et al. 2005; Huard et al. 2006), and L1521F (Bourke et al. 2006), which 
show very different properties in molecular line observations despite their similar internal 
luminosities. IRAM 04191 (Andre et al. 1999; Belloche et al. 2002) and L1521F (Crapsi 



- 2 - 



et al. 2004) show evidence for infall whereas L1014 does not (Crapsi et al. 2005). IRAM 
04191 is associated with a wcll-conimatcd outflow (Andre et al. 1999); the other two are 
not, although at least L1014 and possibly L1521F feature weak, compact outflows (Bourke 
et al. 2005; Bourke et al. 2006). 

The discovery of VeLLOs with Spitzer has put into question the picture of low-mass 

star formation as a continuous process of constant mass accretion at the standard rate of 
2 X IQ-^ Mq yr-^ (Shu, Adams, & Lizano 1987) through a single evolutionary sequence, 
the well-established class system progressing from Class to III (Myers & Lada 1993, Andre 
et al. 1993). This standard accretion rate predicts a much higher luminosity than observed 
for VeLLOs; VeLLOs must feature some combination of a very low central mass and a very 
low accretion rate (e.g. Dunham et al. 2006). If the accretion continues at the current low 
rate to the very small central mass, it might not make a star. However, the accretion rate is 
not necessarily constant. For instance, FU Orionis (FU Ori) objects undergo outbursts (Bell 
et al. 1995 and references therein). Studies (Vorobyov & Basu 2005 and references therein) 
for the nature of the FU Ori variables suggest accretion bursts from the disk to the central 
star by the thermal instability of the disk. 

Therefore, two potential explanations for the very low luminosities of VeLLOs are 1) 
proto-brown dwarfs, and 2) objects in a quiescent phase of the episodic accretion process. 
The former can be discriminated from the latter with studies of the chemistry since they 
involve vastly different thermal histories, which is crucial to the chemical evolution. The 
thermal history is especially important in interactions between gas and ice; ice evaporation 
and gas freeze-out from and onto grain surfaces, respectively, depend on the dust temperature 
(Lee et al. 2004). Proto-brown dwarfs, with their very low masses, will never experience a 
hot phase, whereas the outbursts of a cycle of episodic accretion, a short time period when 
the majority of the mass is dumped onto the central protostar, involve signiflcant warming 
of the surrounding dust. The quiescent states between outbursts feature much colder dust 
temperatures. As a result, envelopes of proto-brown dwarfs will be similar to starless cores in 
their chemical distributions, while objects in a quiescent state of episodic accretion will show 
different chemical distributions from starless cores or normal, embedded Class O/I objects. 

IRAM 04191 may be undergoing episodic accretion since it features a strong outffow 
which predicts a higher accretion rate by two orders of magnitude than inferred from the 
internal luminosity of the source (Andre et al. 1999; Dunham et al. 2006). Furthermore, the 
N2H+ emission, observed with the Plateau de Bure Interferometer (PDBI) and the IRAM 30 
m telescope, shows a hole in the center of the envelope (Belloche & Andre 2004). In general, 
N2H''" emission tends to peak towards the center of starless cores (Lee et al. 2003), but be 
deflcient from the centers of Class O/I sources due to destruction by CO as it evaporates (Lee 
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et al. 2004). IRAM 04191 shows moderate CO depletion (Crapsi et al. 2004); Belloche & 
Andre (2004) suggest that frcczc-out of N2 in the high-density, inner envelope might result 
in the observed N2H"'" hole. However, if freeze-out of N2 is significant enough to explain 
this hole, there should be significantly more deuteration and depletion of CO than observed 
(Crapsi et al. 2004), similar to that seen in prestellar cores (Lee et al. 2003). 

In this study, we model the chemical evolution in the process of episodic accretion to 
provide a possible explanation of the chemical distributions of CO and N2H+ in gas and to 
predict observable consequences in the CO2 ice feature in VeLLOs such as IRAM 04191 that 
show strong evidence for undergoing such a process. 

2. EPISODIC ACCRETION MODELS 

We use the chemo- dynamical model developed by Lee et al. (2004). This model calcu- 
lates the chemical evolution of a model core evolving from the prestellar stage through the 
embedded protostellar stages. The dynamical evolution is described by combining a sequence 
of Bonnor-Ebcrt spheres (Bonnor 1956, Ebert 1955) with the inside-out collapse model (Shu 
1977), where the accretion rate from the envelope onto the star+disk system is constant. 
The model also includes the First Hydrostatic Core (FHSC) stage, which results from the 
first gravitational collapse of a dense molecular core and lasts until the core temperature 
reaches 2000 K and the dissociation of molecular hydrogen causes the second collapse (Boss 
1993, Masunaga et al. 1998). The radius of the FHSC is about 5 AU, and the disk is not 
yet well developed at this stage. As a result, the accretion luminosity arising from accre- 
tion onto the FHSC is not significant. The evolution of the central luminosity follows that 
adopted by Young & Evans (2005), who incorporated the evolution of an unresolved disk 
and stellar photosphere into the central luminosity. This luminosity is proportional to the 
accretion luminosity arising from envelope accretion following the inside-out collapse model 
(see Young & Evans 2005 for details) . 

In this study, we assume that the accretion from the envelope to the disk is constant, 
but the accretion from the disk to the star is episodic. Thus, the envelope density structure 
is identical to that of the inside-out collapse model, but the internal luminosity is changed 
by the episodic accretion process. The initial mass of the model core, with an effective sound 
speed of 0.27 km s~^ and outer radius of 6200 AU, is 1 Mq (Equation 1 of Young & Evans 
2005). The episodic accretion is assumed to begin after the FHSC phase, which lasts for 
2 X 10^ years. 

The dust temperature profile at each time step is calculated with the one dimensional 
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continuum radiative transfer code DUSTY (Ivezic et al. 1999), using as input the density 
profile from the inside-out collapse model. The input internal radiation field at each time 
step is calculated by the model of Young & Evans (2005) with their standard parameters of 
r//) = 0.75, 1]^: = 0.5, and Tmax = 10 (see Young & Evans 2005 for a detailed explanation of 
these parameters) combined with our prescription for episodic accretion. In this prescription, 
accretion bursts 10 times more massive than normal (M ~ 5 x 10~^ Mq yr~^) occur during 
10^ years every 10^ years. For the 9 x 10^ years between the accretion burst events, the 
accretion from the disk to the star is decreased by a factor of 100. Vorobyov & Basu (2005) 
have showed that 0.01 ~ 0.05 Mq is accreted from the disk to the central source with the 
accretion rate of (1 ~ 10) x 10~^ Mq yr~^ during a burst, and the duration of the intervening 
quiescent accretion phase is (1 ~ 3) x 10^ years for M = (1 — 10) x 10~^ M© yr~^. According 
to their calculations and for the accretion rates in the burst and quiescent phases in this 
model, which are about an order of magnitude smaller than the rates of Vorobyov & Basu 
(2005), we adopt the timescales of 10"^ and 10'^ years in the two phases, respectively, and 
these timescales are consistent with the estimated from HL Tauri by Lin ct al. (1994). 

We calculate the radiative transfer independently at each time step since the dust tem- 
perature reaches a new equilibrium essentially instantly compared to the length of our time 
step. For example, one silicate grain with a radius of 0.1 /im and a temperature of 100 K will 
cool down to 10 K within one second if the heating source suddenly disappears (Draine & 
Anderson 1985). Following Young & Evans (2005), the envelope is heated externally by the 
Black-Draine Interstellar Radiation Field (ISRF) (Evans et al. 2001) attenuated by a visual 
extinction of 0.5 mag. Also following Young & Evans, we used the dust model of Ossenkopf 
& Henning (1994) appropriate for thin ice mantles after 10^ years of coagulation at a gas 
density of 10^ cm"^ (0H5 dust). 

Fig. la shows the evolution of the total internal luminosity in the model, and Fig. lb 
shows the temperature evolution of a few infalling grid points. As seen in this figure, the 
temperature that an infalling gas parcel experiences peaks at the time step of each massive 
accretion and is lowered significantly during the quiescent states. Betwen 4 and 5.5x10^ 
years, the luminosity and temperature increase steeply with time even in the quiescent phase 
because of the transition from the FHSC to Class 0, which occurs from 2x10*^ years until 
5.5 xlO'^ years. This transition timescale is the same as that used in Lee et al. (2004), but 
much longer than the 100 years assumed by Young & Evans (2005). 

The evolution of the spectral energy distribution (SED) around 5 x 10"^ years is presented 
in Fig. 2. An accretion burst occurred at 5.1 x 10^ years, after which a quiescent state 
continues until the next burst 10"^ years later. The accretion luminosity powered in the burst 
phase at 5.1 x 10^ years is about 5 Lq, whereas the accretion luminosities at 5.0 and 5.4 
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xlO^ years are about 0.004 and 0.02 Lq, respectively. The model SED at 5.4 x 10^ years 
fits well the observed SED of IRAM 04191. The total luminosity in the model is about 0.09 
Lq, where the luminosity arising from heating by the ISRF (~ 0.07 Lq) is larger than the 
internal luminosity (~ 0.02 Lq). This result is not very consistent with that of Dunham et 
al. (2006), who found that the internal luminosity must be greater than 0.04 Lq. However, 
the density profile and ISRF in their model are significantly different from ours, making it 
hard to directly compare the two results. 

The density and dust temperature profiles are shown in Fig. 3 at the same time steps 
shown for the SEDs in Fig. 2. The infall radius is 0.014 pc (2800 AU) when the accretion 
burst occurs at 5.1 x 10*^. The density at the infall radius is about 2 x 10^ cm~^. The 
density profile has the power law index of ~ 1.5 and 2.0 inside and outside the infall radius, 
respectively. Except for the accretion burst phase, the dust temperature is below the CO 
evaporation temperature (~ 20 K) at almost all radii. 

3. CHEMICAL EVOLUTION IN A CYCLE OF EPISODIC ACCRETION 

In prestellar cores, which are very evolved starless cores, carbon-bearing molecules such 
as CS and CO are readily frozen-out onto grain surfaces (Bergin & Langer 1997, Aikawa et 
al. 2001, Lee et al. 2003, Lee et al. 2004). In contrast, nitrogen-bearing molecules such as 
N2H+ and NH3 are considered good tracers of (column) density in cold, dense cores since N2, 
the parent molecule of and NH3, has a long formation timescale such that it becomes 

abundant later in the core evolution. Therefore, the comparison between CO and N2H"'" 
abundances is a good indicator of the physical properties of a dense core. 

The low binding energy of N2 on grain surfaces has also been suggested as an explanation 
for the long-lived molecules in the cold regions (Lee et al. 2004). However, Obcrg et al. (2005) 
showed that the binding energy of N2 on grain surfaces is not different from that of CO. The 
high abundance of N2H"'' in the cold, dense regions may be due to the replenishment of N2 
through the interactions between N2H+ and grains since N2H+ recombines with electrons on 
the grain surfaces. The significant depletion of CO and electrons, the principal destroyers 
of N2H+ in the gas phase, is another explanation for the high abundance of N2H+ in these 
cold, dense regions. 

We updated the chemical network used in Lee et al. (2004) to include recent results on 
the N2H+ chemistry, including a new binding energy of N2 identical to that of CO (Oberg 
et al. 2005) and new rates of dissociative recombination of N2H+ with electrons (Geppert et 
al. 2004). The initial abundances of other species are the same as those in Table 3 of Lee et 
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al. (2004). Bare silicate grains were assumed for binding energies of molecules. 

Fig. 4 presents the evolution of the CO abundance in gas and ice at the time steps 10^ 
years before, during, and 3 and 5 x 10^ years after the accretion burst at 5.1 x 10^ years 
since the model SED at 5.4 x 10^ years fits the observations of IRAM 04191 reasonably 
well. 10^ years before the accretion burst, CO is frozen-out onto grain surfaces and depleted 
from the gas at most radii except around the CO evaporation radius, inside which the dust 
temperature is higher than the CO evaporation temperature, developed by the previous 
accretion burst at 4.1 x 10"^ years. During the burst phase, the CO evaporation radius 
jumps to a much larger radius due to the significant increase in accretion luminosity. Once 
the burst phase ends and the core makes the transition back to the quiescent phase, the 
accretion luminosity, and therefore the temperature, decreases greatly. CO is then gradually 
frozen back out of the gas and onto the grain surfaces. 

The evolution of the N2H+ abundance is shown in Fig. 5. Before the accretion burst at 
5.1 X 10^ years, the N2H+ abundance peaks at a rather large radius and decreases inward. 
(The sharp drop of the N2H+ abundance at the edge of the model core results from the 
dissociative recombination of N2H+ with electrons.) At radii of 0.002 pc (log r~ —2.7) and 
0.0008 pc (log r~ —3.1), the bumps of N2H+ are caused by the dips of CO at these radii. 
The N2H"'" abundance decreases at radii smaller than 0.0008 pc even though CO is depleted 
here. N2H+ forms mainly by the reaction between N2 and H^, which is destroyed by CO. 
As mentioned above, the evaporation temperature of N2 is the same as that of CO, i.e., ~20 
K. Therefore, at the radii smaller than 0.0008 pc where the dust temperature is lower than 
~ 20 K, and the density is greater than 10® cm^'^, the significant freeze-out of N2 causes the 
low abundance of N2H+. During the burst phase, the N2H+ abundance drops significantly 
at radii smaller than the CO evaporation radius since CO destroys N2H+. In the quiescent 
states 3, 5, and 8 x 10^ years after the burst, N2H+ slowly increases at small radii as CO is 
frozen-out onto the grain surfaces. 



4. DISCUSSION 

Most VeLLOs were classified as starless core before Spitzer observations, and they 
showed wide variation in chemical and physical properties. For instance, IRAM 04191 
(Andre et al. 1999) and L673-7 (M. Dunham 2007, private communication) have well- 
collimated molecular outfiows, but no molecular outfiow has been detected in L1521F even 
though its Spitzer IRAC images show a compact bipolar scattered emission nebula (Bourke 
et al. 2006). L1014 features a very compact outfiow detected with the Submillimeter Array 
(Bourke 2005). In chemical properties, IRAM 04191 and L1014 are similar with moderate 
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CO depletion but different in their N2H+ emission, which peaks at the center of L1014 but 
shows a hole at the center of IRAM04191. L1521F, which was believed to be a highly evolved 
starless core (Crapsi et al. 2004), shows significant CO depletion and centrally peaked N2H"'" 
emission. Before the Spitzer observations, this variation in chemical properties was believed 
to arise from an as-of-yet incompletely understood evolutionary sequence of starless cores. 
The discovery of the VeLLOs has cast doubt on our understanding of the continuous and 
unique sequence of star formation consisting of the evolution of a starless core into a prestel- 
lar core followed by the formation of a Class object that then evolves through the Class 
system. 

The sources without strong outflows might be in the flnal accretion phase of substellar 
objects or in the FHSC stage (Omukai et al. 2007). On the other hand, VeLLOs with strong, 
well-coUimated outflows such as IRAM 04191 are hard to understand since the outflow sug- 
gests a much higher accretion luminosity than the observed internal luminosity of the source. 
A possible explanation for such objects is episodic accretion. In this picture, VeLLOs are ob- 
jects in quiescent phases of a cycle of episodic mass accretion. The strong outflows observed 
can be generated during accretion bursts that occurred prior to the current quiescent phase. 
These two processes predict vastly different thermal histories, which can be discriminated 
by their different predictions for the chemical evolution of the core. Observations probing 
these different predictions are challenging, however, since high spatial resolution is required 
to resolve the regions most affected by the different thermal histories. In addition, observa- 
tions of ice features require high sensitivities due to the very low luminosities of the central 
source. 

As seen in the previous section, the relation between CO and N2H+ abundances in an 
object featuring episodic accretion is different from that of starless cores or hot corinos, 
where CO and N2H+ have an anti-correlation. In the quiescent phase after an outburst, 
both CO and N2H+ have not had enough time to reach chemical equilibrium. As a result, 
CO can be less depleted and N2H"'" more depleted compared to a core forming a proto-brown 
dwarf with a similar central luminosity to IRAM 04191. 

CO in ice can also be a discriminator since CO ice features from VeLLOs in a quiescent 
phase of episodic accretion will be weaker than those in proto-brown dwarfs. Additionally, 
the CO2 ice features can also be very important tracers of the thermal history of VeLLOs. 
The CO2 15.2 iim. bending mode has been detected toward massive star forming regions 
(Gerakines et al. 1999), low mass star forming cores (Pontoppidan et al. 2007) and even dark 
molecular clouds (Knez et al. 2005). The formation process of the CO2 ice is controversial, 
but it must form on grain surfaces since the abundance of the CO2 gas in dark clouds is too 
low (10~^, Bergin et al. 1995) to explain the observed CO2 ice features. CO2 ices are found 
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both in water- and CO-rich grain surfaces. Therefore, the CO2 ice hkely forms on grain 
surfaces by the reactions between CO frozen-out from gas and OH or O photodissociated 
from water or CO in ice. The photodissociation of water or CO on grain surfaces, especially 
in low mas star forming cores, can occur by the cosmic ray induced UV photons. As a result, 
the timescale of the formation of CO2 on grain surfaces must be long (> 10'' years in the 
density of 10^ cm~^, Pontoppidan et al. 2007.) 

The evaporation temperature of the CO2 ice is about 35 K, higher and lower than the 
CO (20 K) and H2O (100 K) evaporation temperatures, respectively. In outburst phases of 
episodic accretion, the dust temperature can rise above the CO2 evaporation temperature 
at some inner radii. There is also a range of radii where the dust temperature is between 
the CO and CO2 evaporation temperatures, so CO2 remains on grain surfaces while CO 
evaporates. In this range, CO2 will produce the double-peaked feature indicative of pure 
CO2 ice (Pontoppidan et al. 2007). In the following quiescent phase, CO will be gradually 
frozen-out as seen above, and thus, combined with the long formation timescale of CO2, the 
solid CO2 abundance profile developed in the outburst phase will not change greatly. As a 
result, we expect that objects in quiescent stages of episodic accretion will have a weaker 
solid CO2 feature (if CO2 evaporates up to a large radius) but a stronger pure CO2 feature 
(if the range of CO2 evaporation is not large) compared to substellar objects. Fig. 6 shows 
the temperature profile during an outburst with marks for the CO and CO2 ice mixture. 
However, if the VeLLOs are proto-brown dwarfs, which never raise the dust temperature 
above the CO and CO2 evaporation temperatures in the majority of their envelopes, then 
the sohd CO2 feature will be strong as seen toward background sources without a strong pure 
CO2 component (Knez et al. 2005). Therefore, the CO2 ice feature as well as the abundance 
distribution of CO and N2H+ in gas can be good tracers of the formation processes of 
VeLLOs. 

This work was supported by the faculty research fund of Sejong University in 2007. We 
are very grateful to Mike Dunham and Neal Evans for valuable comments. 

REFERENCES 

Aikawa, Y., Ohashi, N., Inutsuka, S.-L, Herbst, E. & Takakuwa, S. 2001, ApJ, 552, 639 

Andre, P., Ward-Thompson, D., & Barsony, M. 1993, ApJ, 406, 122 

Andre, P., Motte, F., & Bacmann, A. 1999. ApJ, 513, L57 

Bergin, E.A., Langer, W.D., & Goldsmith, P.F. 1995, ApJ, 441, 222 



- 9 - 



Bell, K.R., Lin, D.N.C., Hartmann, L.W., & Kenyon, S.J. 1995, ApJ, 444, 376 
Belloche, A., Andre, P., Despois, D., & Blinder, S. 2002, A&A, 393, 927 
Belloche, A., & Andre, P. 2004, A&A, L35 
Bergin, E.A. & Langer, W.D. 1997, ApJ, 486, 316 
Bonnor, W.B. 1956, MNRAS, 116, 351 
Boss, A.P. 1993, ApJ, 410, 157 

Bourke, T.L., Crapsi, A., Myers, P.C, Evans, N.J. II, Wilner, D.J., Huard, T.L., J0rgensen, 
J.K., & Young, C.H. 2005, ApJ, 633, L129 

Bourke, T.L., et al. 2006, ApJ, 649, L37 

Crapsi, A., Caselli, P., Walmsley, CM., Tafalla, M., Lee, C.W., Bourke, T.L., & Myers, P.C. 
2004, A&A, 420, 957 

Crapsi, A., et al. 2005, A&A, 439, 1023 

Di Francesco, J., Evans. N.J., II, Caselli, P., Myers, P.C, Shirley, Y., Aikawa, Y., & Tafalla, 
M. 2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil (Tucson: 
Univ. Arizona Press) 

Draine, B.T. & Anderson, N. 1985, 292, 494 

Dunham, M.M., et al. 2006, ApJ, 651, 945 

Ebert, R. 1955, Z. Astrophys., 37, 217 

Evans, N.J.II, Rawlings, J.M.C, Shirley, Y., & Mundy, L.G. 2001, ApJ, 557, 193 

Geppert, W. D., Thomas, R., Semaniak, J., Ehlerding, A., Millar, T., Osterdahl, P., af 
Ugglas, M., Djuric, N., Paal, A., & Larsson, M. 2004, ApJ, 609, 459 

Gerakines, P.A., Whittet, D.C.B., Ehrenfreund, P., Boogert, A.C.A., Tielens, A.G.G.M., 
Schutte, W.A., Chiar, J.E., v an Dishoeck, E.F., Prusti, T., Helmich, P.P., de Graauw, 
T. 1999, ApJ, 522,357 

Huard, T.L., et al. 2006, ApJ, 640, 391 

Ivezic, Z., Nenkova, M., & Elitzur, M. 1999, |astro-ph/9910475 | 

Knez, C, Boogert, A.C.A., Pontoppidan, K.M., Kessler-Silacci, J., van Dishoeck, E.F., 
Evans, N.J. II, Augereau, J.-C, Blake, G.A., & Lahuis, F. 2005, ApJ, 635, L145 

Lee, J.-E., Evans, N.J.,II, Shirley, Y.L., & Tatematsu, K. 2003, ApJ, 583, 789 

Lee, J.-E., Bergin, E.A., & Evans, N.J.II 2004, ApJ, 617, 360 

Lin, D.N.C, Hayashi, M., Bell, K.R., & Ohashi, N. 1994, ApJ, 435, 821 



-10- 



Masunaga, H., Miyama, S.M., & Inutsuka, S.-I. 1998, ApJ, 495, 346 

Oberg K.I., van Broekhuizen, F., Eraser, H.J., Bisschop, S.E., van Dishoeck, E.F., & Schlem- 
mer, S. 2005, ApJ, 621, 33 

Myers, P.C. & Ladd, E.F., 1993, ApJ, 413, L47 

Omukai,K. 2007, PASJ, 59, 589 

Ossenkopf V. & Henning Th., 1994, A&A, 291, 943 

Pontoppidan, K.M, Boogert, A.C.A., Eraser, H.J., van Dishoeck, E.F., Blake, G.A., Lahuis, 
F., Oberg, K.I., Evans, N.J.,II, & Salyk, C. 2007, ApJ, in press 

Shu, F.H. 1977, ApJ, 214, 488 

Shu, F.H., Adams, E.G., & Lizano, S. 1987, ARA&A, 25, 23 
Vorobyov, E.I. & Basu, S. 2005, ApJL, 633, 137 
Young, C.H., et al 2004, ApJS, 154, 396 
Young, C.H., & Evans, N.J. 2005, ApJ, 627, 293 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 11 - 



2 



- (a) 



^ 1 



^ 



-1 



1 \ r 



2x10^ 



1 1 r 



1 1 r 



4x10^ 6x10^ 
time (yrs) 



8x104 



Fig. l.a. — Evolution of the internal luminosity in the episodic accretion model. Outburts 
(massive accretion events) occur every 10^ years for 10^ years after the FHSC phase, which 
lasts until 2 x 10^ years. 
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Fig. l.b. — Evolution of dust temperatures in infalling gas parcels. Grids 200, 300, and 400 
are initially located at the radii of 0.0023, 0.0055, and 0.0124 pc, respectively. 
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Fig. 2. — Evolution of SED in the episodic accretion model. The dashed, thick solid, and 
thin solid lines represent the SEDs 10^ years before, 3 x 10^ years after, and during the 
outburst at 5.1 x 10^ years, respectively. 
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Fig. 3. — Evolution of density and temperature profiles. Time steps for different line styles 
are the same as those in Fig. 2. The density structure in the envelope does not change 
much for 4 x 10^ years, but the dust temperature varies significantly over the episodic 
accretion event. During the outburst, the dust temperature jumps above the CO evaporation 
temperature at radn smaller than about 0.006 pc (~ 1300 AU). During the quiescent states, 
the dust temperature is below the CO evaporation temperature (~ 20 K). 
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Fig. 4. — Evolution of the abundance profile of CO in gas and ice. Time steps for different 
line styles are the same as those in Fig. 2 except for one more thick solid line 5 x 10^ years 
after the outburst event. Before the outburst (dashed line), the CO abundance profiles in 
gas and ice have a peak and a dip, respectively, at ~ 0.002 pc, which are vestiges of the prior 
outburst at 4.1 x 10^ years. During the outburst (thin solid line), CO evaporates, and thus, 
the CO abundance in gas reaches the maximum at radii smaller than 0.006 pc. For the CO 
ice, the opposite happens, i.e., the abundance drops sharply at 0.006 pc. The arrow in each 
panel indicates the evolution of the abundance profile with time after the massive accretion 
event. That is, CO is gradually frozen-out onto grain surfaces to become depleted in gas and 
abundant in ice with time. 
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Fig. 5. — Evolution of the abundance profile of N2H"'". The time steps for different line 
types are the same as those in previous figures, but the thick solid lines show the abundance 
evolution 3, 5, and 8 xlO^ years after the outburst as indicated with the arrow. The bumps 
shown in the dashed line (before the outburst) are vestiges of the previous outburst that 
occurred at 4.1 x 10^ years. During the outburst (thin solid line), N2H''" is destroyed by 
CO inside the CO evaporation radius, so its abundance drops sharply at radii smaller than 
the CO evaporation radius. In the quiescent phase (thick sohd fines), the N2H+ abundance 
gradually increases as CO is frozen out onto grain surfaces. 
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Fig. 6. — Distribution of CO and CO2 ices depending on the temperature distribution 
developed by the outburst at 5.1 x 10^ years. 



